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ABSTRACT 

-/ 

This  paper  is  concerned  with  nonlinear  eigenvalue  problems  of  boundary 
value  problems  for  ordinary  differential  equation  posed  on  an  infinite 
interval.  It  is  shown  that  under  certain  analyticity  assumptions  -  a  domain 
in  the  complex  plain  can  be  identified,  in  which  all  eigenvalues  are 
isolated.  An  intriguing  way  to  solve  such  problems  is  to  cut  the  infinite 
interval  at  a  finite  but  large  enough  point  and  to  impose  additional,  so 
called  asymptotic  boundary  conditions  at  this  far  end.  The  obtained  eigenvalue 
problem  for  the  two  point  boundary  value  problem  on  this  finite  but  large 
interval  can  be  solved  by  any  appropriate  code.  In  this  paper  suitable 
asymptotic  boundary  conditions  are  devised  and  the  order  of  convergence,  as 
the  length  of  the  interval,  on  which  these  approximating  problems  are  posed, 
converges  to  infinity,  is  investigated.  Exponential  convergence  is  shown  for 
well  posed  approximating  problems. 
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SIGNIFICANCE  AND  EXPLANATION 


This  paper  deals  with  the  analysis  and  approximation  theory  of  nonlinear 
eigenvalue  problems  for  boundary  value  problems  of  ordinary  differential 
equation  posed  on  an  infinite  interval.  These  problems  have  the  following 
forms.  We  have  a  linear  homogeneous  system  of  ordinary  differential  equation 
depending  nonlinearly  on  a  complex  parameter  X.  This  system  is  defined  on 
the  interval  tt^,®]  where  tQ  is  some  finite  point.  Moreover  we  have 
linear  homogenous  boundary  conditions  at  tg  and  a  continuity  requirement  at 
infinity.  We  look  for  values  of  the  parameter  X  so  that  a  function  y  £  0 
exists  which  fulfills  the  system  of  differential  equation  over  the  whole 
interval  [ tQ  ,  «>) ,  the  boundary  conditions  at  tQ  and  the  continuity 
requirement  at  infinity. 

Problems  of  this  kind  frequently  occur  in  fluid  mechanics  when  the 
stability  of  a  laminar  boundary  layer  is  investigated  (Orr-Sommerfeld  problem) 
and  in  quantum  physics. 

In  this  paper  properties  of  the  eigenvalues  X  and  the  eigenfunctions 
y  are  investigated. 

An  intriguing  way  to  solve  such  problems  is  to  cut  the  infinite  interval 
at  a  finite  but  large  enough  point  T  and  to  impose  additional  -  so  called 
asymptotic  -  boundary  conditions  at  T,  which  substitute  the  continuity 
requirement  at  infinity.  The  resulting  finite  eigenvalue  problem  can  be 
solved  by  any  appropriate  code,  for  example  by  collocation  methods. 

Suitable  asymptotic  boundary  conditions  are  derived  and  the  convergence 
of  eigenvalues  and  eigenfunctions  as  the  length  of  the  interval,  on  which  the 
finite  approximating  problems  are  posed,  converges  to  infinity,  is 
investigated.  Exponential  convergence  is  shown  for  the  most  important  cases. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  authors  of  this  report. 


NONLINEAR  EIGENVALUE  PROBLEMS  ON  INFINITE  INTERVALS 


Peter  A.  Markowich  and  Richard  Weiss* 


1 .  INTRODUCTION. 

This  paper  is  concerned  with  nonlinear  eiqenvalue  problems  of  the  following  form 

y '  =  taA(t,X)y,  1  <  t  <  »,  a  >  -1  (1.1) 

B( X)y( 1)  =  0  (1.2) 

y  e  C (  [1,o=]  )  :o  y  e  C([1,«>))  and  lim  y(t)  exists  (1.3) 

t  ■>«> 


where  y  is  an  n-vector  and  A(t,X)  is  an  n  x  n  matrix.  Equation  (1.1)  has  a 
singularity  of  the  second  kind  of  rank  a  +  1  in  t  =  ®. 

A  solution  of  (1.1),  (1.2),  (1.3)  is  given  by  a  pair  (p,y),  p  e  C  such  that  y  fe  0 
fulfills  (1.1),  (1.2)  with  X  =  p  and  (1.3).  Eigenvalue  problems  on  infinite  intervals 
occur  frequently  in  quantum  mechanics  and  in  fluid  mechanics,  when  the  stability  of  laminar 
flows  over  infinite  media  is  investigated  (see  Ng  and  Reid  (1980)). 

de  Hoog  and  Weiss  (1980a)  and  Markowich  (1980a)  treated  linear  eigenvalue  problems  on 
infinite  intervals,  i.e.  A(t,X)  =  Aq ( t )  +  XA^t)  and  @  C([1,“])  and  B(  X)  5  B. 

It  was  shown  that  all  eigenvalues  X  of  this  linear  eigenvalue  problem,  for  which  the 
matrix  A(°°,  X)  =  Ag(«o)  +  XA^(«)  has  no  eiqenvalue  on  the  imaginary  axis,  are  isolated,  if 
not  all  X  8  C  are  eigenvalues.  Moreover  if  A^")  =  0  there  is  an  infinite  sequence  of 
eigenvalues  X^  which  fulfill  I X^ I  ♦  ®.  de  Hoog  and  Weiss  (1980a)  also  proved  that  the 
spectral  subspces  are  finite  dimensional. 

The  first  goal  of  this  paper  is  to  show  the  generalization  of  the  isolatedness- 
statement  to  nonlinear  eigenvalue  problems  assuming  that  B(X),  A(t,X)  are  analytic  in 


axis.  The  analyticity  is  supposed  to  hold  for  all  t  @  [1,«]  and  A(t,X)  is  jointly 


continuous  in  x  ft.  The  number  of  rows  of  the  matrix  B(X)  is  assumed  to  equal 

r  ,  which  is  the  sum  of  algebraic  multiplicities  of  eigenvalues  of  A(»,  X)  with  negative 
real  part  for  X  e  SI. 

The  second  goal  of  this  paper  is  to  investigate  the  approximating  eigenvalue  problems 
*4  “  t<XA(t,X)xT,  1  <  t  <  T,  T  »  1  (1.3) 

B(X)x  (1)  =  0  (1.4) 

T 

S(X)xt(T)  *  0  (1.5) 

where  S(  X)  is  a  suitably  chosen  matrix  with  r+  »  n  -  r  rows. 

The  main  question  arising  here  is  to  find  out  which  matrices  S(X)  lead  to 
convergence  of  the  eigenvalues  and  eigenfunctions  of  these  approximating  problems  to  the 
eigenvalues  and  eigenfunctions  of  (1.1),  (1.2),  (1.3)  as  T  ♦  •».  A  class  of  matrices 
S( X)  which  implies  exponential  convergence  will  be  identified.  The  convergence  results 
are  the  generalization  of  the  results  obtained  by  Markowich  (1980a)  for  linear  eigenvalue 
problems.  As  Markowich  (1980a)  pointed  out  there  is  not  always  (even  in  the  case  of  a 
linear  eigenvalue  problem)  an  obvious  way  to  choose  the  suitable  S  which  is  independent 
of  X.  However  there  is  an  intrinsic  way  (see  Keller  (1976))  to  set  up  an  ’asymptotic’ 
boundary  condition  S  depending  on  X.  Therefore  these  'finite'  eigenvalue  problems  are, 
even  in  the  case  of  a  linear  'infinite*  problem,  nonlinear. 

This  paper  is  organized  as  follows,  in  Chapter  2  nonlinear  finite  dimensional 
eigenvalue  problems  are  discussed.  Chapter  3  is  concerned  with  the  case  when  A  is 
independent  of  t,  in  Chapter  4  this  restriction  is  dropped  and  Chapter  5  is  concerned 
with  examples  illustrating  the  theory. 
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2 .  FINITE  DIMENSIONAL  NONLINEAR  EIGENVALUE  PROBLEMS. 

Let  A( X)  be  a  k  «  k  matrix,  holcmcrphic  In  some  domain  n  c  C.  A  value  y  e  C 
for  which  the  linear  equation 

A(y>£  -  0,  i  *  C  (2.1) 

possesses  a  solution,  is  called  an  eigenvalue  and  £  is  a  corresponding  eigenvector.  Let 

detA( X)  denote  the  determinant  of  A(X).  Since  (2.1)  holds  iff  detA(y)  «  0  it  follows 

from  the  identity  theorem  of  holoraorphic  functions  that  either  all  X  6  55  are  eigenvalues 

or  every  compact  subset  of  ft  contains  at  most  finitely  many  eigenvalues. 

Let  e  e  ( 0 , Eg ]  be  a  real  parameter  and  B(X,e)  be  a  k  x  k  matrix,  holomorphic  in 

£1  for  all  e  6  (0,Cg]  and  fulfilling 

lim  sup  IB(X,£)I  *  0  for  all  A  compact,  A  C  il  (2.2) 

£♦0  XGA 

whore  1*1  denotes  some  matrix  norm.  Now  consider  the  perturbed  nonlinear  eigenvalue 
problem 

C(X,£)f  =  (A(  X)  *  B(  X,  e)  )  5  =0.  5  *  o  (2.3) 

Since 

lim  detC(X,E>  «*  detA(X)  (2.4) 

e+0 


-<e  may  employ  standard  perturbation  results  for  zeros  of  holomorphic  functions.  Therefore, 
let  y  be  a  root  of  order  s  of  detA( X)  m  0,  8  be  a  neighbourhood  of  y  and 

b(  c)  »  sup  |detA(X)  -  detC(X,e)l  (2.5) 

xes 

Then  we  get 
Theorem  2.1. 

1  s 

(i)  when  e  is  sufficiently  small  there  are  precisely  s  eigenvalues  u  ,...,y  of 

e  e 

(2.3)  near  y  (counting  multiplicities)  and  they  satisfy 


|UE  - 


<  const,  b (  ep 


. ,  s 


(2.6) 
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N(A(y))  denotes  the  nullspace.  A  proof  can  be  found  in  G.  Vainniko  (1976),  par.  4. 

In  the  case  of  a  linear  eigenvalue  problem  A(A)  E  A  -  XI  we  get  a  stronger 
perturbation  result  for  eigenvectors  if  the  algebraic  and  geometric  multiplicity  of  the 
eigenvalue  u  is  equal  to  one.  Therefore  we  define: 

Definition  2.1.  The  eigenvalue  u  of  (2.1)  is  called  simple  if  u  is  a  zero  of  order 
one  of  detA( X)  =  0. 

It  is  easily  seen  that  p  is  a  simple  eigenvalue  of  (2.1)  iff  det(A(u)  +  xA' ( y) )  =  0  has 
a  zero  of  order  one  at  T  =  0  which  holds  iff  t  =  0  is  an  eigenvalue  of  geometric  and 
algebraic  multiplicity  1  of  the  generalized  linear  eigenvalue  problems 

(A(  \i)  +  TA’(p))£  =  0  .  (2.10) 

We  get 

Theorem  2.3.  Let  u  be  a  simple  eigenvalue  of  (2.1).  Then: 

(i)  for  e  sufficiently  small  there  is  a  unique  eigenvalue  u  of  (2.3)  and  it 

c 

satisfies 

I -  ul  <  const  |detA(u)  -  detc(u,e)|  (2.11) 

(ii)  for  every  u  there  is  exactly  one  eigenvector  £  (with  norm  one)  of  (2.3) 

e  e 

which  satisfies 

•  -  £11  <  const  |detA(u)  -  detc(u,  e)l  •  (2.12) 
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Proof.  The  equations 


* 


C(X,e)£  =0,  15  8  =  1  (2.13) 

e  e 

where  I  •  D  indicates  the  euclidean  norm  in  C  are  a  nonlinear  system  of  equations  for 

(u  r£  )•  The  Frechet  derivative  of  the  unperturbed  problem  (  e=  0)  at  (y,£)  is  given 
£  £ 

by  the  matrix 


(2.14) 


which  is  nonsingular  oecause  y  is  simple.  (i)  and  (ii)  follow  in  a  straight  forward  way 
by  applying  the  techniques  of  Keller  (1975)  and  Vainniko  (1976)  par.  4.  Here  6  shrinks 
to  the  point  y. 

We  conclude  this  section  with  a  result  on  holomorphic  families  of  projections, 
k  k  * 

Theorem  2.4.  Let  P( X)  :  C  C  be  a  family  of  projections  holomorphic  for  X  e  ft.  Then 
04  for  any  pair  (X^Xj)  e  ft  x  ft  there  ia  a  nonsingular  n  x  n  matrix  QfX^.Xj) 
such  that 

P(X,)  =  Q(X1,X2)‘1P(X2)Q(X1.X2)  *  (2.15) 

k  *  k 

(ii)  P(  X.,  )C  is  isomorphic  to  P(  X^)C  for  all  X^Xj 

(iii)  Let  r  =  rank  P(X).  Then  there  is  a  k  x  k  matrix  of  rank  r,  holomorphic  in 
ft,  whose  columns  span  P(X)C  . 

V  * 

Proof.  (i)  follows  from  Kato  (1966)  and  (ii),  (iii)  follow  easily  from  (i). 
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3.  NONLINEAR  CONSTANT  - 

COEFFICIENT  EIGENVALUE  PROBLEMS. 

We  consider 

y*  =  t°A(  A)y,  1  <  t  <  •»,  ci  >  -1 

(3.1) 

B( A)y ( 1 )  -  0 

(3.2) 

y  6  C(  [1,-1  ) 

(3.3) 

where  A(A)  *  0  is  an  n  x  n  matrix. 

The  analysis  for  these  problems  will  outline  the  approach  for  the  more  complicated 
case,  when  A  is  also  a  function  of  the  independent  variable  t.  We  assume  that 
A(*),  B(  •)  are  holomorphic  in  some  domain  ♦  in  the  complex  plain  and  that  there  is  a 
domain  8  C  ♦,  so  that  A(\)  has  no  eigenvalue  v(A)  with  vanishing  real  part. 

Therefore  A(  A)  for  all  A  e  has  affixed  number  of  eigenvalues  with  a  negative  real 
part,  which  we  call  r_  and  a  fixed  number  of  eigenvalues  with  a  positive  real  part,  which 
we  call  r+  (r+  +  r_  »  n).  Now  we  take  a  compact  subset  A  C  0.  Then  there  are  two 
closed  rectifiable  curves  r  ,  T_,  completely  in  the  right  resp.  left  half  plain,  so  that 
for  all  A  e  A  all  eigenvalues  of  A( A)  are  enclosed  by  either  r+  or  T  . 

Now  let 


P(A)  =  tV  f  -  MX))  'dz,  rankP  (A)  “  r  (3.4) 

+  *JTT1  +  + 

P  (A)  “  err  /  (z  -  A(  A) )  'dz,  rankP  (A)  =  r  .  (3.5) 

be  the  total  projections  onto  the  direct  sum  of  invariant  subspaces  associated  with 
eigenvalues  of  A ( X )  with  positive  resp.  negative  real  part. 

Prom  Kato  (1966),  Chapter  2  we  conclude  that  P+,  P_  are  holomorphic  in  A',  the 
interior  of  A. 

The  general  solution  of  the  problem  (3.1),  (3.3)  is 

od-1 

y(t,A)  =  exp  (  ^  -  a(X)>JAH,  c  e  cn  (3.6) 
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Theorem  2.4,  (iii),  implies  that  there  is  a  n  x  r  matrix  V( X)  of  full  rank  and 
holomorphic  in  A*  which  spans  P  (X)Cn.  Therefore  we  can  rewrite  (3.6)  and  get, 
inserting  into  the  boundary  condition  (3.2) 

F(X)n  2  B(X)exp  y)  v(  X)  n  -  0,  T)  e  C  ~  (3.7) 

assuming  that  B(X)  is  a  r  x  n  matrix.  Every  pair  (ij,n),  n  *  0  which  solves  (3.7) 
determines  a  solution  of  the  eigenvalue  problem  (3.1),  (3.2),  (3.3)  by 

cd-1 

y(t,u)  =■  exp  (--+-  A(  u)  )V(  |i)  n  .  (3.8) 

Our  assumptions  guarantee  that  F(X)  2  B(X)exp  V(X)  is  holomorphic  in  A“,  so  we 

get  from  Chapter  2 

Theorem  3.1.  Let  B ( X )  be  a  r  x  n  matrix.  Then  either  all  X  6  0  are  eigenvalues  of 

(3.1),  (3.2),  (3.3)  or  every  compact  subset  of  fl  contains  at  most  a  finite  number  of 

eigenvalues.  If  u  is  an  eigenvalue  of  (3.1),  (3.2),  (3.3),  the  dimension  of  the 
nullspace  is  between  1  and  r_. 

Now  we  want  to  approximate  the  eigenvalue  problem  (3.1),  (3.2),  (3.3)  by  finite 
interval  problems: 

x'  =  t°A( X)x  ,  1  <  t  <  T,  T  »  1  (3.9) 

T  T 

B(X)X  (1)  =  0  (3.10) 

T 

S( X)x  (T)  =0  (3.11) 

T 

where  S( X)  is  a  r+  x  n  matrix  whose  choice  will  be  discussed  later. 

We  write  the  general  solution  of  (3.9)  as 

a+1  gtl 

x  ( t )  =  exp  f-  —  -  A(  X)  )V(  X)n  +  exp  f- - - -  A(X)]w(X)n  (3.12) 

T  va  +  1  -  v  a  +  1  '  + 

where  the  columns  of  the  n  x  r+  matrix  W( X),  which  is  holomophic  in  A°,  span  P+(X)cn. 
Evaluation  at  the  boundaries  and  using  (3.10),  f3.11)  gives  the  n  x  n  block  system 
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f 


e  i \,T) 


B (  X)exp  X) 

+  1 

or*- 1 


nGt-  1 


B( X)exp 


a  +  1 


A (X) ;W( X) 


S(X)exp  (•= - -A(X)V(  X)  S(  X)W{  X) 

a  +  1 


1  r 


j 


O  (3.13) 


We  conclude  that 

detF(X,T)  *  detF(X)  •  detS( A)W( X)  +  c(X,T>  (3.14) 

holds,  where 

_  *  +  1  1 

— T  i' 

|c(  X,T)  |  <  const.  Uexp  '■ - -  A(  X)  ]w(  X)  tt  •  U  S  (  X)exp  ( - r  A(  \)  V( X)  «  (3. IS) 

v  a  +  i  '  v  &  +  i 

Let  v  (X)  be  the  largest  negative  real  part  of  the  eigenvalues  of  A ( X )  and  let  v  ( X) 
be  the  smallest  positive  real  part  of  the  eigenvalues.  Then  (3.15)  reduces  to 


_a+ 1 


I  c  (  X,T)  I  <  const  (  X)  exp  ,(v  (X)  -  v  (  X)  ) 


(3.16) 


a  +  1  * 

where  const(X)  is  bounded  when  X  varies  in  a  compact  set.  Now  we  prove  the  convergence 
theorem. 

Theorem  3.2:  Let  the  r+  x  n  matrix  S( X)  be  holomorphic  for  X  e  ft  and  assume  that 

detS(X)W(X)  *  0  for  X  8  0.  Let  u  6  ft  be  an  eigenvalue  of  (3.1),  (3.2),  (3.3)  of 

order  s,  i.e.  detF(X)  has  a  zero  at  X  =  u  of  order  s.  Then  there  are  exactly  s 
1  s 

eigenvalues  uT<...<UT  (counting  multiplicities  of  the  zeros  of  detF(X.T))  for  T 
sufficiently  large  in  a  sufficiently  small  neighbourhood  of  u  and 


i  .  t 

max  Ivi  —  u I  <  const,  exp  v  (u)  -  v  (  u >  +  e)  — — — — ' 
i=  1  ( 1  )s  '  sin  i 


(3.17) 


I  Uj,  -  ill  <  const,  exp  ((v_(u)  -  v+(y)  +  e)  -  ^  1 


( 3. IP) 


1  ?  i 

where  ^  =  —  >  li^  and  e  >  0  is  sufficiently  small.  Let  xT  be  an  eigenf unct ic 


i=  1 


belonging  to  one  of  the  U^,'s.  Then 


Tn+1 

inf  lx  -  y  n  <  const,  exp  ((u  (  u)  -  y  (  v)  +  c)  — - - -  1  (3.  IP) 

yeN  11,1  ”  S<“  11 

u 


where  N  denotes  the  nullspace  of  (3.1),  (3.2),  (3.3)  for  X  =  u* 
U 


Proof .  All  statements  follow  immediately  by  regarding  (3.13)  as  perturbed  eigenvalue 
problem  of 


B(A)exp  [-tii-)v(X) 
va  +  1  ' 


SU)W(X) 


(3.20) 


(which  is  equivalent  to  (3.7))  and  by  applying  the  Theorems  2.1  and  2.2. 

Now  we  discuss  a  possible  choice  of  S(X).  Let  the  rows  of  the  r  x  n-matrix  S  (X) 

+  P 

T 

span  the  range  of  (P^(A))  (the  superscript  denotes  transposition).  Then  the  asymptotic 
boundary  condition 

S  ( A)x  (T)  =  0  (3.21) 

P  T 

fulfills  the  assumptions  of  Theorem  3.2.  Moreover 


P  +  (  X)exp  A(  X)  )v(  X)  5  0,  X  e  SJ 


(3.22) 


holds.  Therefore,  using  the  boundary  condition  (3.21),  the  matrix  in  the  (2.1)  position  in 
(3.13)  vanishes  for  all  T  >  1  and  the  approximate  problems  (3.9),  (3.10),  (3.11) 
reproduce  the  eigenvalues  and  eigenfunctions  of  the  problem  (3.1),  (3.2),  (3.3)  exactly. 


4.  GENERAL  NONLINEAR  EIGENVALUE  PROBLEMS  ON  INFINITE  INTERVALS. 

We  consider  the  following  problem 

y’  =  taA(t,X)y,  1  <  t  <  «,  a  >  -1  (4.1) 

B( A)y( 1 )  =  0  (4.2) 

y  S  C(  [1,»]  )  (4.3) 

where  A(t,X)  is  an  n  x  n  matrix  holomorphic  for  A  in  some  domain  and  every  fixed 
t  e  (1/®]  and  continuous  in  [ 1 , <*>)  x  0.  Also  B  is  holomorphic  in  We  assume  that 

there  is  a  domain  ft  C  i ,  so  that  the  matrix  A{  X)  has  no  eigenvalue  v( A)  on  the 
imaginary  axis  for  X  6  ft.  As  in  Chapter  3  we  take  any  compact  subset  A  ft  and 


construct  the  projections  P  (X),  P  ( X) 


(4.4) 


P .  C  X )  =  -r— 7  f  (z  -  A(«,  X))  'dz,  rankP  (  X)  =  r 


2ni 


(X)  =  —7  /  (z  -  A(  «»,  X )  )  'dz,  rankP  (  X)  S  r 
-  2ni  i  - 


The  contours  r  ,  T  are  chosen  as  in  Chapter  3.  We  set 


...  , .  cA(<»,  X)  .  a+1  < 

$(t,X)  =■  exp  . - —  t  j 

^  a  +  1  ‘ 


(4.5) 


(4.6) 


and  define  the  operator  :  C([6,H)  •>  C((6,<»!)  for  6  >  1 

(H  g)(t)  =  f  P+(X)(^'1(s,X)sag(s)ds  + 

CD 

t  (4.7) 

+  <t>(  t ,  X  3  /  P  (  X>  <t>~ 1  (s,  X,sag(  s)ds 

<5 

so  that  H^g  8  C([6,«l)  is  a  particular  solution  of  the  problem 

y'  =  taA(  <»,  X)y  +  tag(t),  t  >  5,  g  e  C!  [«,»)  )  .  (4.R) 

An  analysis  of  can  be  found  in  de  Hoog  and  Weiss  ( 1980a, b).  P+ ,  P  are  holomorphic 

in  A®  for  every  (fixed)  t  e  and  continuous  for  t  e  (1,®).  Then  it  is  an  easy 

exercise  to  show  that  (H  g(«,X))(t)  is  holomorphic  in  A®  for  every  fixed  t  e  [6,°°]. 

From  de  Hoog  and  Weiss  ( 1980a, b)  we  conclude  that 

AH  II  i  <  C(X)  (4.9) 

X  l  6, 

where  C( X)  is  independent  of  5. 

Now  we  show  that  C(X)  remains  bounded  when  X  varies  in  compact  subsets  K  C  A®. 
From  (4.7)  we  derive 

00  t 

IH  ».  <  max  (/  IF  !t,s,  X)  lids  +  /  IF  (t,s,XHds)  (4.10) 

1  te  [{,<»)  t  s' 


where 
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(a)  F+(t,s,X)  -  s%(t,X)P+(XH_t(s,X),  XCA 

<b)  F_(t,s,X)  -  a%(t,X)P_(X)*',(9,X),  XCA 
holds.  We  transform  A(»,  X)  to  its  Jordan  canonical  form  J(»,  X): 

a(»,x)  -  e(X)j<«,x>e~1(x>,  xen 

and  assume  that  has  the  block  structure 

J(«,  X)  -  diag( J+(  X) , J~(  X)) 

OB  OD 

where  the  r+  x  r+  matrix  J*(X)  contains  only  eigenvalues  with  positive  real  part 
r_  x  r_  matrix  J^(X)  contains  only  eigenvalues  with  negative  real  part  for  all 
Defining  the  diagonal  projection 

D+  -  diagdr  ,0),  D_  -  diag(0,Ir  ) 

we  get 


P+  (  X )  -  E(X)D+E_1(X),  P_<X)  -  E( X)D_E*\ X),  let 


(4.11) 

(4.12) 

(4.13) 
and  the 
X  e  0. 

(4.14) 

(4.15) 


and  obtain 


(a) 


F  (t.s.X)  -  s^fX) 


J+(  X) 


exp  (ttt  (t 


or+1  a+1 


))  0 


e"\x) 


(b) 


F  (t,s,X)  -  s^tX) 


J  (  X) 


exp(7TT  (t 


cr* 1  ct+ 1 


)) 


E_1(  X) 


Obviously,  F  (t,s,*)»  F  (t,s,»)  are  holomorphic  in  A0  for  fixed  s,t. 
Each  entry  of  F  ,  F  is  a  sum  of  the  form 


(4.16) 


f  (t,s,X)  *  sa  V  a  (X)exp  f— - -  (t 

i  i 


v,  (  X )  .  .  .  .  j 

oc+1  a+1.->,.o+1  a+1.  i 

-  s  )  )(t  -  s  ) 


(4.17) 


where  v,  (X)  are  the  eigenvalues  of  J  ( X)  and  v. ( X)  are  the  eigenvalues  of  J  (X). 
i.  w  i  °° 
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The  integers  fulfills  0  <  <  (r+  -  1)  resp.  0  <  j  <  (r  -  1).  a^(X)  is  a  sum 

of  products  of  elements  of  E(X)  and  E  '(X). 

Now  we  take  a  compact  subset  K  c  A°.  It  follows  from  Kato  (1966)  that  E(X',  E  1 ( X ) 

can  be  chosen  holomorphically  in  K  -  (z z  }  where  the  z.  are  points  at  which 

in  i 

eigenvalues  v*(X)  change  algebraic  or  geometric  multiplicities.  Also  J(X)  is 

holomorphic  in  K  -  {z  ,  ...,z  }.  Without  loss  of  generality  we  assume  that  none  of  the 
1  n 

Zj^'s  lies  on  the  boundary  3K.  If  that  happens  we  can  choose  a  larger  compact  set  K  D  K, 

so  that  . . .  }  n  3K  =  $.  Since  the  entries  f  (t,s,X)  are  holomorphic  in  K,  they 

in  ± 

take  their  maximum  at  the  boundary  3K.  The  coefficients  a,  (U  and  the  v*(A)  are 

i  l 

continuous  on  3K  and  therefore 


max  | f  ( t ,s,  X) |  < 


max  |a  .  (  X)  |  max  (exp( — -  (t**1  -  s^1. 


..  ct+1  cttl  i  . 
It  -  s  | 


Xe3K 
j=1 ( 1 )r 


where  c.  =  min  Re  u.  (X),  c  =  max  Re  v  ( X)  hold.  Therefore  we  get 
1  Xe3K  X  1  Xe3K  1 


(a)  max  HF  (t,s,X)R  <  c  (K)sa  max  (expf - -  (t  -  s  )l(s  -  t0^  )x]  (4.19) 

+  I  d+l  1 

XeK  0  <i <r 


(b)  max  IF  (t,s,  X)  «  <  const.  sa  max  (exp( - -  (t3*1  -  sa+1))(t°f1  -  s3*1)*}  (4.ig) 

XeK  0<i<r  3  + 


will  0  <  c  =  min  c  ,  0  >  c  =  max  c.  and  0  C  i  <  n  hold. 

+  i=1(1)r+  1  *  i= 1 ( 1 ) r  _  1 

Using  (4.10)  and  the  estimates  derived  in  Markowich  ( 1980b)  we  get 


max  IK  I  ,  <  C(K) 
XeK  X  !5'“’ 


where  C(K)  is  independent  of  6. 


We  rewrite  (4.1)  as 


setting 


We  get  from  (4.21) 


y'  “  taA(~,\)y  +  ta(A(t,X)  -  A(«,X))y 


G(t,  X)  =  A(t,X)  -  A(»,X) 

r 

y(t)  -  <{>(t ,  X)V(  X)  n  +  (HjG(  •/  X)y)  (t) ,  gee 


where  the  columns  of  the  (n  x  r  )-matrix  V(X),  which  are  holomorphic  in  A*,  span 
range  of  P_(X).  The  assumptions  on  A(t,X)  guarantee  that  there  is  a  5  >  1,  6  « 

such  that 

IG(  • ,  X)  1 .  <  — for  all  XeK 

2C(K) 

where  K  is  any  compact  subset  of  A“  and  C(K)  is  as  of  (4.20).  Therefore 


max  |H  G( 

xex  A 


X>,(«,»)  4  2 


This  implies  that  I  -  H^G(*,X)  :  C(  t  <5,  >  ♦  C(!6<'»))  is  nonsingular  for  all  X  e  K 

-1  r- 
y  =  (I  -  HXG(.,X))  ♦(•,X)V(X)n,  gee 


holds,  y  is  defined  for  t  e  (6,“]  and  all  X  e  K.  The  series  expansion  for  the 
n  x  r  matrix 

<|i_  ( t ,  X )  =  ((I  -  HXG(  •,  X)  )_1  $(  •,  X)V(  X)  )  (t) 

is  given  by 

00 

^(•,X)=  l  (H  G{  .,X))%(  .,X)V(  x)  e  C(  [  5,»1 ),  xex 
i=o  A 

( k )  ° 

The  partial  sums  ifi  (t,A)  of  this  series  are  holomorphic  in  X  e  K  for  all  fixed 

t  e  1 6 » and  because  of  (4.25)  we  get 


(4.21 ) 

(4.22) 

(4.23) 

the 

6(K), 

(4.24) 


(4.25) 

and 

(4.26) 


4.27) 

4.28) 
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(4.29) 


li|/k)vt)l  <  (  l  max  1 0(  •,  X)V(  X)  I 


i=0 


X6K 


(6,-) 


$(t,X)V(X)  is  holomorphic  in  A°  0  K  and  therefore  the  partial  suras  are  uniformly 
bounded  on  K  and  so  (t,X)  is  holomorphic  in  X  e  £  for  all  fixed  t  @  By 

continuation  (t,X)  is  holomorphic  in  X  S  K  for  all  fixed  t  S 

Inserting  into  the  boundary  condition  (4.2)  gives  the  finite  dimensional  eigenvalue 
problem 

r 

F(X)n  =  b(  X)  i(i_<  i,  X)  n  =  o,  nec"  (4.30) 

where  B( X)  is  assumed  to  a  r  x  n  matrix.  Given  now  any  compact  subset  8  C  8  we 

0 

choose  A,  K  such  that  A*  D  K,  K  D  0.  So  F(  X)  is  holomorphic  in  9  and  Theorem  3*  1 
holds  for  the  problem  (4.1),  (4.2),  (4.3).  Therefore,  excluding  the  trivial  case,  all 
eigenvalues  in  ft  are  isolated  and  the  dimension  of  the  nullspace  is  between  1  and  r_. 
Now  we  prove  the  asymptotic  estimate  for  ^  (t,A): 


max  l\\>  (t,X)B  <  const.  exp((  max  Re  v.  (  X)  +  e)  — — r) 

xee  "  Aeae  “ 

i*  1  ( 1 )  r 


cr+1 


(4.31) 


where  ( X )  are  the  eigenvalues  of  A<»,  X)  with  real  part  less  than  zero  and  e  >  0  is 
sufficiently  small  so  that  the  exponent  has  negative  sign. 

If  there  is  one  (or  more)  of  the  singularities  {z  ,...,z  }  of  E( X)  on  the  boundary 
36,  we  take  a  larger  set  9^  such  that  9  C  9^ »  9^  C  K  and  {z^,...,z^}  O  361  =  $ . 

Then  we  derive  as  in  (4. IP) 


t0^1 

max  «4>(t,X)8  <  const.  exp((  max  Re  v.  (X)  +  e)  •  — 

xee  xeag  1  a 

i=  1 ! 1 )  r 


(4.32) 


The  necessity  to  add  e  >  0  in  the  exponent  comes  from  the  possibility  that  0  might  have 
to  be  changed  to  8  described  above  and  from  the  possible  occurence  of  powers  of 

It®*'  -  s®*1 | .  A  sufficiently  small  change  and  the  continuity  of  the  eigenvalues  assure 


that  c  is  sufficiently  small. 
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max  I (H  G(-,X)$(*,X))(t)l  < 
X08 


<  max  IG(  •, X)  1 .  .  .  (/  max  IF  (t,  s,  X)  i  max  l$(s , X)  Ids  + 

xee  t  xee  +  xee 


t 

+  /  max  IF  (t,s,X)l  max  l#(s,X)Ids) 

&  X08  '  X0  8 

Now  (4.19),  (4.32)  can  be  used  for  the  estimation  of  the  righthand  side  of  (4.39). 
to  be  substituted  for  K  in  the  definition  of  c*,  c^.  Since  e  >  0  the  estimate 
in  Markowich  (1980b),  Chapter  1,  Theorem  2.3,  can  be  used  and 


max  l(H  G(  •,  X)4><  *,X)(t)  I  4 
X80  x 

4  const,  max  IG(*,X)I.,  .exp((  max  Re  \>  ( X)  +  e) 
X08  1  X0  38 

i- 1 ( 1 )r_ 


o+1 

h-r» 


follows.  Repeated  use  of  (4.34)  and  (4.28)  gives  (4.31) 

As  in  Chapter  3  we  investigate  the  approximating  ’finite’  eigenvalue  problems: 
x^  -  taA(t,X)xT,  1  4  t  4  T,  T  »  1 
B(X)x  (1)  -  0 
S(X)xt(T)  -  0 

where  S(X)  is  a  suitably  chosen  r+  x  n  matrix  whose  entries  are  holomorphic  in 
Rewriting  (4e35)  as 

x’  =  t°*A( »,  X)x  +  taG(t,  X)x  ,  1  4  t  4  T 

T  T  T 

where  G(t,X)  is  defined  as  of  (4.22),  we  set 

xT  =  <Mt,x)v( X)n_  +  *(t,x)*'1(T,x)w(X)n+  +  (»x  TG(  •,x)xT)(t) 

where  the  columns  of  the  n  x  r+  matrix  W(X)»  which  can  be  chosen  holomorphic  In 
span  the  range  of  P+(X)  and  ^  :  C(I6,T1)  *►  C([6,T])  is  defined  as 

HX,T9  =  HXgT 

for  g  e  C( (6,T) ),  1  4  6  4  T  and 


(4.33) 


6  has 
given 

(4.34) 


(4.35) 

(4.36) 

(4.37) 

0. 

(4.38) 

(4.39) 
A”, 

(4.40) 
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qT(fc>  ■ 


g(t),  5  <  t  <  T 
g(T) ,  t  >  T 


(4.41) 


has  been  set. 

Given  a  fixed  eigenvalue  X  »  p  e  fl  of  (4.1),  (4.2),  (4.3),  we  take  a  compact  subset 
K  C  /\«  with  u  6  K  and  conclude  from  (4.20), 

MX  »H  <  max  1H,I..  ,  <  C(K)  (4.42) 

xex  X'T  t4'T1  xex  x  t*'"1 

and  therefore  there  Is  a  fixed  6  *  6(K)  >  1  such  that 


max  «HA(TG(.,X).{5t]  <- 


(4.43) 


holds  so  that  (I  -  H  G(*.X))  1  exists  for  all  X  e  K  as  an  operator  on  C((fi,T]).  We 

Af  T 

get  from  (4.39)j 

x_  -  (I-H,  _G(  .,X))"\( -,X)V(  X)h  +  (I-H.  _G(  *,X))"%(  * ,  X)  ij,-1  (T ,  X)WC  X)h  (4.44) 

T  A/T  —  A/T  + 

on  (S,“l.  The  analyticity  of 


(a) 

(b) 


T*_(t,X)  -  ((X  -  H ^  TG(*,X))-1*(*,X)V(X))(t) 
>|>+(t,X)  -  ((I  -  H ^  TG(-,X))"%(«,X)i)i"1(T,X)W(X))(t) 


(4.45) 

(4.45) 


in  X  for  t  e  { <5 , T]  follows  as  the  analyticity  of  f  (t,X). 

The  n  x  r  matrix  ^  resp.  the  n  x  r+  matrix  iji+  fulfill  the  equations 


(a) 

(b) 


Tt_(*,X)  -  <HX  TG(-,X)Ti|»_<*,X»)(-)  -  <{>(  - ,  X  )V(  X ) 
i|i+<  • ,  X)  -  (H^  TG(  *,X)T\|i+(  •,  X)H  •)  -  <t>(  •,  X)  <t>~  7  (T,  X)W(  X) 


(4.46) 

(4.46) 


Similar  to  de  Hooq  and  Weiss  (1980a)  we  derive  some  properties  of  From  (4.27) 

and  (4.46)  (a)  we  get  for  TZ_  =  i|i_  -  if 


Z  =  H  G  (  • ,  X ) _Z  +  (H  _G {  - ,  X )  if  -  H  G(  • ,  X)  if  ) 

X  “  A/T  T  ”  A/T  ”  A  - 


(4.47) 
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and  therefore  we  get  by  regarding  G(  •,  X)  iji_(  *.  X)  e  C([6,T]) 


• ,  X)  -  4>_(  • ,  X ) 


(I  - 


H,  G(*,X))-1(H,  _  -  H  >G(  • ,  X)  4>  (-.X)  e  C(  t  6,T]  ) 
X»  T  A#  T  A  “ 


Obviously  for  g  e  C([6.T))  and  t  e  [6<T] 

T  . 

((H,  ~  "  H  )g)(t)  =«  <J>(  t ,  X)  /  P  JXH  (s,X)s“(g(T)  -  g(s))ds  = 

A#T  A  oo  + 

T 

=  4lt,\)f  (T,X)P+(X)$(T,X>  /  P+(X)$“  (s,X)sa(g(T)  -  g(e))ds  - 

00 

»  4.(t(X)*-1(T,X)W(X)y(g,T) 


where  y(g,T)  e  C  +  holds  and  therefore 

<(H  ~  “  H  )G(«,X)t|i  (vX))lt)  -  *(t,X)*'1(T,X)W(X)r 

A»T  A  “  1 

is  fulfilled.  r  is  a  r  X  r  matrix.  From  (4.47)  and  (4.48)  we  derive 
T  + 

T> *  <!/_(•.  X)  +  T^+(  *,  x>rT 

Therefore  the  matrix  [\|i  (t,X),  Ti|i+(t,A)]  has  rank  n  for  all  t  e  [6,T]  and  is 
fundamental  matrix  of  (4.35). 

Instead  of  (4.44)  we  can  write  the  general  solution  of  (4.35)  as 

xT  «  iji_(t,x)n_  +  Ti(i+(t, X) n+ 

For  the  following  we  need  an  estimate  for  TI('+<*.X)'  From  (4.45)  (b)  we  derive 

max  lTi|»+(*.X)  -  *(  *,  X)^"1  (T,  X)W(  X)  l[6  TJ  < 

XG  0 

<  const,  max  IG(  •,  X)  $(  »,  X)  $  1  (T,  X)W(  X)  I .  , 


Using  similar  analyticity  arguments  as  above  it  is  easy  to  check  that  the  right  hand 
of  (4.52)  can  be  estimated  by 

w(T,k)  =  max  (max  «G(  t ,  X)  lexp (— - ^  (t0*1  -  T3*1)] 

tets.Ti  xee  a 


(4.48) 


(4.49) 


(4.50) 


(4.51) 


(4.52) 


side 


(4.53) 
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J 


where  k  5  k( 6)  ”  min  Re  v*(A)  and  e  >  0  sufficiently  small. 

xeae  1 

i-1(1)r+ 

Obviously  lim  w(T,k)  »  0  and  we  get  after  continuation  to  [1,T] 

T+co 


lim  .  X>  -  ♦(.,U^'(T,l)liUH  =  0  (4.54) 

T*»  1,1 

uniformly  for  A  6  9. 

Now  we  evaluate  (4.51)  at  the  boundaries  t  =■  1,T  and  use  the  boundary  conditions 
(4.36),  (4.37)  getting  the  n-dimensional  nonlinear  eigenvalue  problem 


n 


B(  A)  ij)_(  1 ,  A) 
S{  A)  i|/_ (T,  A) 


B(A)T*+(t,A) 
S(A)t«|.+  (T,A) 


n 


Interpreting  (4.55)  as  perturbed  eigenvalue  problem  of 


(4.55) 


B(A)*_(1,A)  0 

0  S(A)W(A)  +  o(T,  A) 


(4.56) 


we  get  with  F(  A)  ■  B(  A)  \|i_(  1 ,  A) : 


detF(  A,T)  «  detF(  A)  (detS(  A ) W (  A)  +  o(T,X))  + 
+  0(1  i|p+(1,A)«  •  DS(  A)  \|i_  (T,  X)  I ) 


(4.57) 


where 

I o (T, A)  I  +  0  as  T  ♦  “  (4.58) 

uniformly  for  A  S  8. 

Assuming  that  S(X)W(X)  is  nonsingular  we  get  by  cancelling  detS(A)W(A)  +  o(T,X) 
and  by  applying  the  perturbation  arguments  of  Chapter  2: 

Theorem  4.1.  Let  the  r+  x  n  matrix  S(  A)  be  holomorphic  for  )  e  SI  and  assume  that 
detS(A)W(A)  #0  in  0.  Let  u  8  (5  be  an  eigenvalue  of  (4.1),  (4.2),  (4.3)  of  order 
s,  i.e.  detF(A)  =*  detB(A)\|i  ( 1 ,  A)  has  a  zero  of  order  s  at  A  *>  y.  Then  there  are 
exactly  s  eigenvalues  u^,,...,u^  (counting  multiplicities  of  the  zeros  of  detF(A,T)) 
for  T  sufficiently  large  in  a  sufficiently  3mall  neighbourhood  S  of  p  and 


i 
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(4.59) 


r 


„<*i 


max  I  wi  -  ul  <  const. (w(T, x(S  ))]8exp((v  (y)  +  e)  — p - — ) 

1  M  '  S(a  ’ 


where  v  ( y)  ■  max  Re  v,(y)  and  w(T,x)  is  defined  in  (4.53)  and  e  >  0  is 
i«1(1)r_ 

sufficiently  small. 

„0rH 


(4.60) 


I  UT  -  w!  4  const.  w(T,  x(S^)  )exp((  v_(  y)  +  e>  ?  —  —  ) 

*  1  r  i  i 

Hg  m  ~  1  Wj  has  *5een  set.  Let  xT  be  an  eigenfunction  belonging  to  one  of  the  n^’s. 


Then 


.or*- 1 


inf  lx  -  yl,  ,  <  const.  fw(T,x(I  ))1  exp((v  (  y)  +  e)  — — — )  (4.61) 

T  l  ■  $  *  J  M  ~  810  ♦  I  ) 

VJ 


holds  where  N  denotes  the  nullspace  of  (4.1),  (4.2),  (4.3)  for  X  •  y. 
y 

These  convergence  results  are  the  extension  of  the  convergence  results  for  linear 
eigenvalue  problems  given  in  Markovich  (1980).  The  orders  of  convergence  obtained  there 
hold  without  any  change  for  nonlinear  problems. 

A  possible  choice  for  S(X)  is  given  by  (3.21),  i.e.  the  rows  of  the  holomorphic 

T 

r+  x  n  matrix  SIX)  •  S^( X)  span  the  range  of  (P+(X))  (the  superscript  T  denotes 
transposition).  This  choice  reproduces  eigenvalues  and  eigenvectors  exactly  in  the  case 
that  A  does  not  depend  on  t.  However,  in  the  general  case  this  does  not  hold  anymore 
but  in  some  Important  cases  the  asymptotic  boundary  condition  Sp(X)xT<T)  «  0  implies  a 
faster  order  of  convergence  than  given  in  Theorem  4.1.  Assume  that  A(t,X>  decays 
algebraically  or  exponentially: 

Aft,  X )  -  A(X)  +  0(tVa(t))  for  t  +  »  (4.62) 

uniformly  in  compact  subset  K  C  JJ  where  -y  e  R  and  a(t)  >0  is  a  real  function  such 
that  t^e  ***"'  *o  as  t  ♦  «.  Then  since  S  ( X)  t(T,  X)V(  X)  =  0  we  get  from  (4,46a) 

IS  (  X)  ( T ,  X )  I  -  IS  (  X)(HG<  XH  (*,X))(T)I  < 

P  -  p  X 

y  -a(T)  .  1 

<  const.  T'e  exp[(  max  Re  v.(X)  +  c)  - - \ 

i"1  ( 1  )r  1  a 
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This  follows  from  the  estimates  given  in  Markowich  (1980b)  applied  to  (4.27).  In  this  case 

the  right  hand^side  of  the  estimate  (4.58),  (4.61)  given  in  Theorem  4.1  can  be  multiplied 

y  -a  (T)  a  v 

by  (T'e  )  and  the  right  hand  side  of  (4.60)  can  be  multiplied  by  T'e  ’ . 

Now  we  consider  the  case  of  simple  eigenvalues  of  (4.1),  (4.2),  (4.3).  Since  we  only 

have  defined  simple  eigenvalues  for  nonlinear  finite-dimensional  eigenvalue  problems  we 

give 

Definition  4. 1 :  An  eigenvalue  y  e  ft  of  (4.1),  (4.2),  (4.3)  is  called  simple  if  the 
corresponding  nullspace  is  one  dimensional,  say  it  is  spanned  by  the  normed  vector  y,  and 
if  the  problem 


v*  -  t<JA(t,y)v  -  t°^  (t,y)y(t) 

(4.63) 

B(  y)v(  1 )  +  B  (u)y(1)  -  0 

(4.64) 

v  e  c(  [i,»]  > 

(4.65) 

has  no  solution.  Now  we  show 

Theorem  4.2.  The  eigenvalue  y  of  (4.1),  (4.2),  (4.3)  is  simple  iff 

y  is  a  first  order 

zero  of  detF( X)  ”0. 

Proof:  Since  y  is  an  eigenvector  corresponding  to  the  eigenvalue  y  e  ft 

y(t,y)  » 

r 

for  some  (  e  C  holds.  Obviously 

yX<t'u>  “  fx 

is  a  particular  solution  of  (4.63),  (4.65)  therefore  the  general  solution  of  (4.63)  is 

r 

v(t)  *  \|Mt,vi)0  +  yx(t,y),  gee* 

Inserting  into  (4.64)  gives 

B(ll)l|>_(1#u)B  *  “*  ( B  (  y  )  4>_(1,u)  +  B  x  (  y )  ijj  (  1 ,  y )  )  £ 

or 

F(u)B  =  -Fx(uU  • 

This  equation  is  unsolvable  (for  B)  iff  the  generalized  linear  eigenvalue  problem 
(F(y)  +  <F x ( y ) )  *  0  has  <  »  0  is  an  eigenvalue  with  geometric  and  alqebraic 

multiplicity  1.  This  holds  iff  det  F(X)  *  0  has  a  first  order  zero  at  X  *  y. 


-20- 


can  be 


Now  we  show  that  approximations  for  an  eigenvalue-eigenvector  pair  (X,y) 
computed  -  in  the  case  that  \  is  simple  -  as  solutions  to  nonlinear  'finite'  two  point 
boundary  value  problems.  Lentini  and  Keller  (1980)  did  computation  pursuing  this  way. 

We  set,  assuming  that  yen  is  simple 


Yn+1  -  V.  *  ‘  <V, 


(4.66) 


and  get  from  (4.1),  (4.2),  (4.3) 


z' 


n 


t“f(t,z),  1  <  t  <  » 


(4.67) 


*,m\ 


iV1*' 


I  V"  -  1 

i*1 


»  b(z(  1 ) )  *  0 


(4.68) 


z  e  C(  [!,•)) 


(4.69) 


The  condition  j  y  ( 1 )  «  1  shall  sort  one  eigenfunction  y  t  0  out  of  the  one 
i=1  1 

dimensional  eigenspace. 

(4.671,  (4.68),  (4.69)  is  a  singular  two  point  boundary  value  problem  as  described  by 
de  Hoog  and  Weiss  ( 1980a, b),  Markowich  ( 1980a, b,c)  and  Lentini  and  Keller  (1980).  Since 

all  eigenvalues  in  (1  are  isolated  and  since  u  is  simple  the  solution 
T 

z  =  (y,,...,y  ,w)  of  (4.67),  (4.68),  (4.69)  is  locally  unique.  Now  we  will  show  that  z 
1  n 

is  isolated,  i.e.  the  linearized  problem,  has  only  the  zero-solution.  We  get  for  the 

T 

linearized  problem  with  u  =  <u u  ,u  , ) 

1  n  n+1 


u' 


A(  t ,  y) 
0 


A^(t,u)y(t) 

0 


u 


(4.70) 
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r 


7 


B(u> 


2y(  1) 


B  i  u)Y< 1) 


u(1|  =  0 


(4.71 ) 


u  e  C(  [1,»]  ) 


Setting  v=  (u  ,...,u  )  we  derive 
1  n 


(a) 
<b) 
(a) 
( b  > 


v‘  -  t°X(t,  iilv  =  u  4taA,(t,u)\ 
n+ 1  A 


n+1 


-  const. 


B!u>v(1)  +  u  ( 1 )B  ( y)y( 1 )  =  0 
n+1  A 


(4.72) 

(4.73) 

(4.73) 

(4.74) 

(4.74) 

(4.75) 


y<1)  v( 1 )  =  0 

v  e  C<  [1,»]  )  • 

Because  of  Definition  4.1  the  problem  (4.73),  (4.74),  (4.75)  has  no  solution  unless 
u^+1  *  0.  If  un+1  *  0  then  v  has  to  be  an  eigenfunction  of  (4.1),  (4.2),  (4.3), 
therefore  v  =  cy  for  some  constant  c.  (4.741(b)  gives  c  =  0,  such  that  u  =  0 
follows  as  the  unique  solution  of  (4.70),  (4.71),  (4.72).  Therefore,  we  conclude  from 
Markowich  (1980b)  that  the  infinite  problem  (4.67),  (4.68),  (4.69)  can  be  approximated  by 
finite  interval  problems  of  the  form 

(4.76) 

0  (4.77) 

(4.78) 


t  f(t,w  ),  1  <  t  <  T 


b(wT( 1 ) ) 


S  ( w  ( T ) )  =  0 
T 


where  w^  =  (w^, . . .  ,w^,w^  1  )T  and  (w^,,...,w”)  is  the  approximation  to  the  eigenvector 
y  and  w™  1  is  th<=  approximation  to  A.  The  choice  of  s  :  Rn+'  ♦  R  +  is  explained  in 
Markowich  (1980b).  The  analysis  given  here  shows  that  we  can  take 


S(w  (T) ) 
T 


S 

P 


(«' 


n+1 


(T)  ) 


(4.79) 


The  superscript  T  denotes  transposition.  Markowich  (1980b)  showed  that  the  solution 

wT  is  locally  (around  z)  unique  for  T  sufficiently  large  and  that 

Iz  -  w  I  <  const.  1S(z (T) )  1  =  const.  US  (p)y(T)l  (4.R0) 

T  (1 ,T1  p 
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holds  and  therefore  we  get  the  order  of  convergence  given  in  Theorem  (4.1)  with  s  -  1, 

because  the  boundary  condition  (4.78)  is  equivalent  to  (3.21). 

If  (4.63)  holds  the  order  of  convergence  is 

I U_  -  u I  <  const.  T^exp((  max  Re  v  (X)  +  e)  -  a(T)) 

^  i-1(1)r_ 

and  the  same  holds  for  the  normed  eigenvectors. 


5.  CASE  STUDIES. 

The  first  problem  we  treat  is  the  so  called  radial  Schrodinger  equation  of  the  Kepler- 
problem  (see  JUrgens-Rellich  (1976),  Chapter  3,  par.  9)  which  is  given  by 

-u"  +  {(1  +  Ulr  2  -  2cr  1  }u  =  -Xu,  1  4  r  <  »  (5.1) 


where  ten,  c  e  R  holds. 
0 

The  transformation 

takes  (5.1)  into  the  system 


A(r,  X) 


4  r  <  » 


(5.2) 


(5.3) 


such  that 

0  1 

A(-,x)  =  (5 

X  0 

holds. 

The  eigenvalue  problem  (5.3)  is  linear  but  we  will  construct  the  nonlinear  (in  X) 
asymptotic  boundary  condition  S^(X). 

The  Jordan  form  J(»,X)  of  A(»,X)  is 


'  /x  o  1  0  1  on] 

J(»,X)  =  _  ,  r  a  1,  r_  -  1,  D+  -  ,  D  =  (5.5) 

o  -A  J+  “  +  l  o  o  J  •  [o  1  J 
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and  therefore  the  set  £1  '{IS  C|Re(/x  )  *  0}  is  given  by 


£1  =  C  -  {X|Re  X  <  0}  .  (5.6) 

With  appropriate  boundary  conditions  of  the  form 

a,y 1 ( 1 >  +  a2y2(1>  =  0  <5. 7) 

we  conclude  from  Chapter  4  that  every  eigenvalue  X  6  £1  of  (5.3),  (5.7)  is  isolated  and 
that  the  dimension  of  the  nullspace  equals  1. 

A  complete  analysis  of  the  problem  is  given  in  Jurgens  and  Rellich  (1976).  They  show 
that  there  is  an  infinite  sequence  of  eigenvalues  X*°' 

X(n)  =  c2(i  +  1  +  n)"2  V  n  e  Nq  (5.8) 

and  the  eigenfunctions  y^n'  are  given  by 


(n),  ,  -  ( 

y  (r)  =  exp[-y 


.  (n) 


v  1+1  ,  , 

:)r  Pn(r) 


exp(-/X 


(n) 


)ri+n+1(1 


0  ( r  ' 1 )  ) 


(5.9) 


because  pn(r)  is  a  polynomial  in  r  of  degree  n.  They  assumed  that  a^  =  sin  a, 
a^  »  cos  a  with  a  e  (0,°°). 

A  straightforward  calculation  gives  for  X  e  Q 


1 

1 

1  1 

[  /X  -/X  1 

- 1 

2/X 

2X 

2  JTx 

E<  X)  = 

,  E  '(X)  - 

,  (  X )  = 

lx  > 

1 

1 

+ 

/x  1 

.  -i7x 

2  X 

.22 

and  therefore,  since  P^(X)  *  E(X)D+E  1 ( X) 

S  (X)  =  t/x,1)  (5.11) 

P 

holds. 


The  approximating  problems  have  the  form 

x*  =  A(r,X)x  ,  1  <  r  <  R,  R  >>  1 

R  R 

fai '*2,XR< =  0 


(5.12) 

(5.13) 


[ /X , 1 ] x  (R)  =  0~  /X  xl ( R)  +  xZ(R)  =  0,  x  =  (x1,x2)T 
R  R  R  R  R  R 


(5.14) 
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1 


Since  A(r,X) 

X  0 

analysis  given  in  chapter  4  show9  that 


+  0(r  ')  and  since  w(R,Re/x  )  <  const.  R  1  the  convergence 


i,,(n)  . (n) . 

max(|X  -  X^  | 


,  (n)  (n)  •.  -Jx' n ^  R  i+n-1 

,y  -  XR  '  I  1 .  R]  <CS  R 


(5.15) 


under  the  assumption  that  the  X^n'’s  are  simple.  The  XRn*’s  resp.  xRn^’s  are  the 

eigenvalues  resp.  eigenvectors  of  (5.12),  (5.13),  (5.14)  which  approximate 

,(n)  (n) 

X  resp.  y 

(5.15)  follows  by  using  \(r  (r,X^n'l  »  (y^n'(r),  y^n'(r)]^.  Therefore  the  estimate 

-  or 

(4.81),  where  e  appears  in  the  exponent,  can  be  improved  using  (4.57).  From  plarkowich 


(1980a)  we  conclude  that  every  boundary  condition  SxR(R)  =  0  where  S  is  independent  of 


X  and  where 


SE ( X ) C  Z  )  *  0  in  a 


(5.16) 


holds,  leads  to  convergence  of  the  order 


(IXtn)  -  4n’l,.y(n)  -  x'n>«(1,Rj)  * 


(5.17) 


<  c1 1y*n'  (R)  lw(R,Re(  /x  ))  <  c  exp(-/xln)  r)r 


(n)  f+n 


Setting  S  «  (s^Sj!,  e  C,  (5.16)  is  fulfilled  iff 

®1  +  S2^  *  °  (5‘ 181 

holds.  For  example  the  natural  boundary  condition 

xR( R)  =0  (5. 19) 

fulfills  (5.18)  (s^  *>  1,  s  «  0)  and  the  order  of  convergence,  given  by  (5.17)  differs 
only  by  one  power  of  R  from  the  order  of  convergence  produced  by  the  ’optimal1  boundary 
condition  (5.14). 
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The  second  problem  we  deal  with  is  the  Orr-Sommerfeld  equation  (see  Ng  and  Reid 


(1980))  which  governs  the  stability  of  a  laminar  boundary  layer  in  a  parallel  flow 
approximation . 


—  (■ 

iRa  ' 


dz 


a2}2*  -  {(U(z)  -  X)  (— r  -  a2)$  -  U"(z)*|  -  0 
dz 


(5.20) 


with  a  >  0,  R  >  0  is  the  Reynolds  number/  U(z)  is  the  velocity  distribution  fulfilling 


U(z) 


1  +  F(z)e-li;Z  ,  U  >  0,  F  e  C2((0,»n 


(5.21) 


such  that  U(<°)  »  1,  U"(<»)  =  0  holds.  4>(z)e^a*X  represents  the  disturbance  stream 

function.  The  boundary  conditions  for  the  Orr-Somraerf eld  problem  are, 

<M0)  =  ( 0 )  *  $(«)  =  *•(»)  «  0  (5.22) 


The  substitution 

y  =  (*,*•, 4", 4>"’)T  (5.23) 

gives  the  linear  eigenvalue  problem 

0  1  0 
,  0  0  1 

y  0  0  0 

f^z)  0  f2(z) 

* - 

A(  z ,  X) 


+  X 


0  0  0  0 

0  0  0  0 

0  0  0  0 

a  0  b  0 


0  4  z  <  od 


(5.24) 


where 

(a) 

(b) 


10  0  0 

y( 0 )  =  0  (5.25) 

a  i  o  o 

y  e  C(  (0,»]  )  (5.26) 

f ,  ( z )  =  -(a  +  iaR(a2U(z)  +  U"(z>)>  (5.27) 

f  2  ( z )  =  2a2  +  iaRU(z)  (5.27) 


(c) 

<d) 

holds. 

V,U>  -  a,v2U)  -  (a2  ♦  iaR(1-X))1/2  ,v3(X)  -  -a,v4(X)  -  -(a2  +  ioR(1-X))1/2 


ia3R 


-iaR 


The  eigenvalues  of  A(<°,  X)  are 


so  that  Re  v.(X),  Re  v_(X)  >  Oi  Re  u,(X).  Re  v.(X)  <  0  for  all 
12  3  4 

X  0  n  ■  C  -  {X|Re  X  ■  1,  I«  X  <  •  ^}.  Xll  eigenvalues  X  e  fi  of  (5.24),  (5.25),  ( 
are  isolated  and  the  nullspacea  are  at  most  two-dimensional. 

The  approximating  problems  have  the  form 

x^  -  A(z,  X)x  ,  0  <  z  <  Z 


0  0 
1  0 


0 

0 


xz(0>  -  0 


S(X)xz(Z)  -  0 


12  3  4 

where  x  »  (x  , x  ,x  ,x  )  holds  and  X  0  fl. 

z  z  z  z  z 

As  for  the  first  example  we  calculate  the  'optimal'  boundary  condition  S(X>  = 


S  (X) 
P 


0  1  v2<X)  av2(X> 

1  0  -(v2<X)  +  av2<X!  +  a2)  -av^(  X)  ( \i^(  X)  +  a) 


Since  A(z,X)  «  A( X)  +  0(z2e  WZ  }  holds  we  get  from  (4.81)  for  simple  eigenvalues 
X  -  u  0  fi 


maxllu  -  Uzl,ly  -  xz'[0fZ,>  < 


<  const  Z3e  J,Z  w(Z,min(a,Re  v2  (  y) )  )exp(max( -a, Re  v^fy)  +  e)z) 


(5.27) 

(5.27) 

(5.28) 

5.26) 

(5.28) 

(5.30) 


S  (X) 

P 

(5.32) 


(5.33) 
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where  y,xz  are  the  norraed  eigenfunctions.  In  the  most  interesting  case 

2  2 

a  <  |Re  ( \ ) !  <  1  the  order  of  convergence  is  Z  exp(-wZ  -  2(a  -  e)Z),  t  >  0 

sufficiently  small,  and  linear  asymptotic  boundary  conditions,  for  example 

x\z)  =  x2(Z)  »  0  (5.34) 

z  z 

achieve  a  slower  order  of  convergence,  namely  exp(-2(a  -  e)Z)  (see  Markowich  (1980a)). 

For  this  problem  the  nonlinear  asymptotic  boundary  conditions  achieve  a  much  faster 
order  of  convergence. 

Numerical  experiments  were  performed  by  Ng  and  Reid  (1980). 
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